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O ■ Abstract 

in 

We consider staircase polygons, counted by perimeter and sums of k-th pow- 
ers of their diagonal lengths, k being a positive integer. We derive limit distribu- 
tions for these parameters in the limit of large perimeter and compare the results to 
Monte-Carlo simulations of self-avoiding polygons. We also analyse staircase poly- 
gons, counted by width and sums of powers of their column heights, and we apply 
our methods to related models of directed walks. 



dedicated to Tony Guttmann on the occasion of his 60th birthday 



1 Introduction 



Self-avoiding polygons are (images of) simple closed curves on the square lattice. We are 
interested in counting such polygons by perimeter, enclosed area, and other quantities, 
where we identify polygons which are equal up to a translation. In statistical physics, self- 
avoiding polygons, counted by perimeter and area, serve as a model of two-dimensional 
vesicles, where the vesicle pressure is conjugate to the polygon area [21 HH]- Despite the 
long history of the problem (see fl§\ HH] for reviews), there is not much exact, let alone 
rigorous knowledge. Recent progress comes from two different directions. It appears that 
stochastic methods, combined with ideas of conformal invariance, describe a number of 
fractal curves appearing in various models of statistical physics, including self-avoiding 
polygons [21]. This will not further be discussed in this article. 

Another line of research consists in analysing solvable subclasses of self-avoiding poly- 
gons. Such classes of polygons have been studied for a long time and by different people, see 
e.g. [7] and references therein for (subclasses of) column-convex polygons. In this article, 
we will concentrate on one such subclass, the class of staircase polygons. The perimeter 
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generating function of staircase polygons has been given by Levine and by Polya [32] , 
the width and area generating function has been given by Klarner and Rivest [23J an d by 
Bender [I], the perimeter and area generating function of this model has been given by 
Brak and Guttmann [Oj, Lin and Tseng j2E], Bousquet-Melou and Viennot 6], Prellberg 
and Brak |3I|, and by Bousquet-Melou [7. , using various solution methods. 

Of particular interest to the Melbourne group of statistical mechanics and combinatorics 
is the critical behaviour of such models, in particular their critical exponents and scaling 
functions. For staircase polygons, counted by perimeter and area, these quantities had 
been obtained in 1995 [33J LHj On the basis of the numerical observation that rooted 
self-avoiding polygons, counted by perimeter and area, have the same critical exponents as 
staircase polygons, Tony Guttmann raised the question [35] whether their scaling functions, 
which describe the vesicle collapse transition, might be the same. One of the implications 
of this conjecture concerns the asymptotic distribution of polygon area in a uniform model, 
where polygons of fixed perimeter occur with equal weight. This distribution was tested 
for by self-avoiding polygon data from exact enumeration and found to be satisfied to a 
high degree of numerical accuracy |37j. A field-theoretic justification for the form of the 
scaling function was then given by Cardy ^U]. Later, other implications of the scaling 
function conjecture relating to the asymptotic distribution of perimeter at criticality and 
to the crossover behaviour were found to be satisfied, within numerical accuracy |%H] , 

One might thus be tempted to conclude that critical self-avoiding polygons and critical 
staircase polygons lie in the same universality class. Then, one should expect that also 
distributions of parameters different from the area should coincide for both models. Inves- 
tigating the implications of this idea is one of the aims of this article. We thus re-analyse 
staircase polygons, which we count by a number of parameters generalising the area. In 
so doing, we adopt the viewpoint of analytic combinatorics [12!, thereby highlighting con- 
nections to statistical physics and to a probabilistic description. For example, staircase 
polygons are in one-to-one correspondence to Dyck paths fTJ H3J , see Figure ^ for an illus- 
tration and below for definitions. In turn, Dyck paths are classical objects of enumerative 




Figure 1: The set of staircase polygons is in one-to-one correspondence with the set of Dyck 
paths ^n[33]- Vertical column heights of a polygon correspond to peak heights of a path. 

combinatorics [I3J , which are asymptotically described by Brownian excursions. This con- 
nection is described at various places in the literature, see e.g. [221 HH EE] Moreover, Dyck 
paths are closely related to models of trees JBII3]> as indicated in Figure [2] In fact, the 



2 



o o 

o 

Figure 2: The set of rooted, ordered trees is in one-to-one correspondence with the set of Dyck 
paths ^2133]. Leave distances to the root of a tree correspond to peak heights of a path. 

generating functions of polygon models, models of (simply generated) trees, and models 
of paths satisfy similar functional equations, typically leading to limit distributions of the 
same type. More generally, it appears that limit distributions for counting parameters are 
irrelevant of many details of the functional equation. This phenomenon of universality 
has recently been investigated using the setup of g-functional equations |H] and extends 
previous results [HI EHj by including counting parameters which generalise the polygon 
area. The results in 01] were obtained by an interplay between methods from statisti- 
cal physics, combinatorics and stochastics. Corresponding counting parameters appeared 
previously within a combinatorial setup, see Duchon 14 s and Nguyen The |rU)[ \M\. The 
derivation of moment recurrences of limit distributions for these parameters is inspired by 
the method of dominant balance as used by Prellberg and Brak [31] in polygon statistical 
physics, and some results for the corresponding stochastic objects appear to be unknown 
in the stochastics literature. 

The outline of this article is as follows. We first derive a functional equation for the 
generating function of staircase polygons, counted by perimeter and sums of k-th powers of 
lengths of diagonals, k being a positive integer. We also analyse a corresponding quantity 
for column heights. We then derive limit distributions for the above counting parameters 
in the limit of large perimeter, resp. width. We will first review the case k = 1, which serves 
as a preparation for the following section, where the case of general k is discussed. We 
then investigate the question of universality of our results by comparing them to Monte- 
Carlo simulations of self-avoiding polygons. The simulations indicate that the considered 
parameters have different laws for self-avoiding polygons and for staircase polygons. The 
methods discussed in this paper can be applied to similar problems concerning models of 
walks and models of trees. This we indicate for various models of discrete walks related to 
Dyck paths, such as meanders and bridges. The analysis extends previous work [3TJ1 I3T] , 
where mainly the case k = 2 is studied. 
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2 Functional equations 



Consider two fully directed walks on the edges of the square lattice (i.e., walks stepping 
only up or right), which both start at the origin and end in the same vertex, but have no 
other vertex in common. The edge set of such a configuration is called a staircase polygon, 
if it is nonempty. Each staircase polygon is a self-avoiding polygon. For a given staircase 
polygon, consider the construction of moving the upper directed walk one unit down and 
one unit to the right. For each walk, remove its first and its last edge. The resulting 
object is a sequence of (horizontal and vertical) edges and staircase polygons, see Figure 
El The unit square yields the empty sequence. Using this construction, it is easy to see 



Figure 3: The set of staircase polygons is in one-to-one correspondence with the set of ordered 
sequences of edges and staircase polygons. A corresponding bijection is geometrically charac- 
terised by shifting the upper walk of a staircase polygon one unit down and one unit to the left, 
and by then removing the first and the last edge of each walk. 

that the set of staircase polygons V is in one-to-one correspondence with the set Q of 
ordered sequences of edges and staircase polygons. Let us denote the corresponding map 
by / : V — > Q. Thus, for a staircase polygon p e V, we have f(p) — (q±, . . . , q n ) E Q, where 
is, for % — 1, . . . , n, either a single edge or a staircase polygon. The image of the unit 
square is the empty sequence n = 0. A variant of this decomposition will be used below in 
order to derive a functional equation for the generating function of staircase polygons. 

The horizontal (vertical) perimeter of a staircase polygon is the number of its horizontal 
(vertical) edges, its (total) perimeter is the sum of horizontal and vertical perimeter. The 
horizontal half-perimeter is also called the width of the polygon, the vertical half-perimeter 
is also called the height of the polygon. The area of a staircase polygon is the number of its 
enclosed squares. The half-perimeter of a staircase polygon is equal to the number of its 
(negative) diagonals plus one, the area of a staircase polygon equals the sum of the lengths 
of its (negative) diagonals, see Figure |U for an illustration. For k a nonnegative integer, we 
will consider k-th diagonal length moments of a staircase polygon p, defined by 



where the summation ranges over the set of diagonals D(p) of p, with 1(d) denoting the 
length of d, i.e., the number of squares crossed by d. We count staircase polygons by half- 
perimeter no(p) + l and by the parameters nk(p), where k = 1, . . . , M. Let u = (m , . . . ,%) 






d£D{p) 
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Figure 4: A staircase polygon of half-perimeter 14 and area 18. It has 8 diagonals of length one 
and 5 diagonals of length two. 



denote formal variables. For a staircase polygon p, define its weight w p (u) by 

W P (U) = U n / P)+ \T (P) ■ ■ ■ ■ ■ u n M MiP \ 

and let G(u) = J2 pe -pW p (u) denote the generating function of V. The bijection described 
above can be used to derive a functional equation for G{u). 

Theorem 2.1. The generating function G(u) of staircase polygons satisfies 

G ( u ) - u 2 oUl -...-u M , 2) 

where the functions Vk(u) are given by 

M (') 

= n«i (A; = 0,1,...,M). (2.3) 

i=k 

Proof. For a sequence q — (q±, . . . q n ) G Q, its preimage p = /~ 1 (q r ) G V can be described in 
the following way, which will prove more convenient for enumeration purposes. Consider 
the sequence (pi, ■ ■ ■ ,p n ) of staircase polygons, where Pi = f° r * = ^,---,n. 

The staircase polygon p is obtained from . . . ,p n ) by concatenation, i.e., by translating 
such that the top right square of Pi and the lower left square of p i+ \ coincide, for 
i = 1, ... ,n — 1. Let us denote by V C "P the subset of staircase polygons such that for 
each p G V we have p = where g is either a single edge or a staircase polygon. The 

above description implies that the set of staircase polygons is in one-to-one correspondence 
to the set of concatenated ordered sequences of staircase polygons from V . 

For a single edge e, the weight of the staircase polygon / _1 (e) is given by Wf-i( e )(u) = 
Uqu\ ■ . . . ■ u 2 M . For a staircase polygon p, note that by the binomial theorem we have 

n k {r\p)) = E = 2 + E (^) + !)* = 2 + E 

d6D(/-i(p)) den(p) i=0 ^ ' 
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This implies for the weights 

M M / I \ / M \ / M M 

«/- W («) = % n «,- (r,w) = «» n «? n ^ = «. n «? n n 

i=0 Z=0 \ fc=0 / \«=0 / \fc=0 l=k 

(M \ ( M \ (M M N 

=< n«' n«'«r w nn«. u "* <p) 

\«=1 / V=0 / \fc=l «=fc 

(M \ M 
JJm; ) W (^) JJ^fc(^)" fe(p) = u\ux • • • • • MmW p (u(«)), 
Z=l / fc=0 

where the functions v^{u) are given by (|2.3|) . For a sequence (pi, . . . ,p n ) of staircase 
polygons, denote their concatenation by c(pi, . . . ,p n )- We have ty c (p 1)P2 )(tt) = Uq 2 u^ 1 ■ . . . ■ 
u~^w Pl (u)w P2 (u) , a monomial in -u, and thus 

w c(pu ..., Pn) (u) = 2(n _ 1} — —^p^) ■■■■■w Pn (u). 

u u 1 ■ . . . ■ u M 

The above bijection yields a functional equation for the generating function G(u). We 
have 



G{U) = ^ Yl W c(pl,-,Pn)( U ) 

n=0 (pi,...,p n )G(P)™ 

oo 

= MqUi • . . . - u M } j 



w pi (u) w Pn {u) 



n=0 (Pl,-,Pn)e(P)" 




UqWi • • • • • mm 



1 — 2uqUi ■ . . . ■ um — G(y{uj) 
This is the functional equation (|2.2j) . □ 



Remark. The case M = is standard, see e.g. jjo]. The case M — 1 has been studied 
by various techniques in [T%1 ITT] [7]. For corresponding models of walks and trees, see 
also the review in ^7j. Analogous counting parameters for the case M > 1 appear, within 
the context of wall polyominoes, in [T3j. This study has been extended in [3U1 |3T] and 
in j^T] . Whereas the focus in |3UJ EI] is on particular examples of directed walks such as 
Dyck paths, |HJ (resp. also discuss the structure of underlying classes of functional 
equations in the case M = 1 (resp. M > 1). Explicit expressions for the generating function 
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G(u) have been obtained for M = and M — 1, see the above references. A functional 
equation for staircase polygons, counted by the above parameters and by width and height 
in addition (with formal variables x and y), can be derived by the same method. It is 
obtained from ()2.2j) by multiplying the rhs with xy and replacing the factor of two in the 
denominator by (x + y) . 

Instead of considering diagonal length moments, one can also analyse k-th column height 
moments of a staircase polygon p, defined by rrikij)) = J2 c ec(p) K c ) k i where the sum is over 
all (vertical) columns c of p, with C (p) denoting the set of (vertical) columns of p, and with 
1(c) denoting the number of squares in column c. Obviously, m (p) is the width, mi(p) the 
area of p. Counting polygons by column height moments, width and height h(p) leads one 
to consider weights 

w p (u,y) = y h ^.ur ip) -...-uZ M{P \ 

y being a formal variable. Let H(u,y) = ^2 p& pW p (u,y) denote the corresponding gen- 
erating function. For a staircase polygon, consider the construction of moving the upper 
walk of the polygon one unit down, and removing its first edge and the last edge of the 
lower walk. Similarly to the case discussed above, this construction can be used to derive 
a functional equation for H(u,y). It is 

H(u, y) = y H(v(u)y)+u .....u M 

1 - (H(v(u),y) + u ■ . . . ■ u M ) 

where the functions Vk(u) are again given by (|2.3|) . 

Consider now the problem of counting staircase polygons by (total) half-perimeter and 
by column height moments. The generating function of the model is given by H(u,uq). 
Since the sum of column heights of a staircase polygon is its area, for M = 1 we have 
H(uq, u±, Mo) = G(uq, ui) with functional equation ()2.2j) . For M > 1, the diagonal moment 
approach for deriving a functional equation cannot be adapted to the present case, since 
in the inflation step, which includes the addition of a vertical layer, the height of a column 
may be increased by more than one unit. It is not possible to keep track of this varying 
height increase without introducing new auxiliary variables in the generating function. 
But then, the method of extracting limit distributions of section 0] cannot be applied in 
its present form. Since we prefer a model where limit distributions are considered with 
respect to large perimeter (and not to large width), we chose to discuss the diagonal length 
moment model. However, this choice is no restriction, since limit distributions of the first 
two models are of the same type, see section |3J Also, an exact analysis of the first few 
moments of the third model suggests the same type of limit distribution, see section El 
below. 



3 Perimeter and area generating function 

In this section, we will review the case M — 1, which corresponds to counting polygons by 
half-perimeter and area. Whereas most results discussed below have appeared at various 
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places in the literature j^U El El EH EH] , we recollect them here informally in order to 
prepare for the following section, where the case of general M is discussed. 

Consider the functional equation (|2.2J) for M — 1 and set uq = t and u\ = q. The gen- 
erating function G(t, q) of staircase polygons, counted by half-perimeter and area, satisfies 
the functional equation 

The function G(t) := G{t, 1) is the generating function for staircase polygons counted by 
half-perimeter. If q = 1, the functional equation (j3.1J) reduces to a quadratic equation, 
whose relevant solution is given by 



G (t) = 1 2t -^O* = — ^ f 2n V n+1 = t 2 + 2t 3 + 5t 4 + 14t 5 + 0(t 6 ), 

n=l ' 

where C n = ( 2 ^) /(n + 1) is the n-th Catalan number. Thus, the number p n of staircase 
polygons of half-perimeter n is given by 

1 /2n — 2\ 1 4 n_1 



The above asymptotic form for p n can alternatively be obtained by singularity analysis of 
the generating function G(t). Such an approach will also be used in the following section. 
The function t > G(t) is analytic for \t\ < t c = 1/4 except at t = t c . The above expression 
yields an expansion about t — t c of the form 

oo 

G(t) = G(t c ) + J2ki(tc-tY l+1)/2 , (3.2) 

2=0 

where G(t c ) = 1/4, / ,o = — 1, io,i = 1, and /o,/ = for / > 1. In order to obtain 
information about the coefficients of G(t), we remind that functions g(x) = (x c — x) -7 
with rational exponent 7 ^ { — 1, —2, . . .} satisfy 

[x n ]g(x) = ^^fc^- 1 (1 + Oin- 1 )) , (3.3) 

where [x n ]f(x) denotes the coefficient of order n in the power series f(x), and T(z) denotes 
the Gamma function [1211121 • Applying this result to the (finite) expansion ()3.2|) yields 
the above asymptotic form for p n . 

We are interested in a probabilistic description of staircase polygons. Let V n C V 
denote the subset of staircase polygons of half-perimeter n. Consider a uniform model, 
where each polygon p e V n occurs with equal probability. We ask for the mean area of 
a polygon of half-perimeter n. To this end, introduce a discrete random variable X n of 
polygon area by 

P(X„ = m) = 

Z^m Pn,m 
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where p n , m denotes the number of polygons of half-perimeter n and area m. Note that the 
numbers p n>m appear in the expansion G(t,q) = J2 n m Pn,mt r ' l q m of the perimeter and area 
generating function G(t,q). The mean area of a staircase polygon of half-perimeter n is 
then given by 

^ \t n ] q^-G(t, q) 



E[X n ] 



9=1 



E m Pn,m [t n ]G(t, l) 

The function appearing in the numerator can be obtained from the functional equation 
(J3.1|) by differentiating w.r. t. q and setting q = 1. This yields 4-G(t, q)\ = t 2 /(l - At), 

q q=l 

and we get 

E[X n ] = £1 = ^ 2 (1 + 0(n-i)) . 

Thus, the mean area of a staircase polygon scales with its perimeter as r?l 2 . For a single 
staircase polygon, the mean length of its diagonals is given by the quotient of area and 
half-perimeter n minus one. Thus, the mean length of a diagonal of a random staircase 
polygon scales with its perimeter as n 1 / 2 . The mean width of a staircase polygon may be 
analysed similarly using the width and height generating function, which is algebraic of 
degree two, see ([2.4)1 . It can be shown that the mean width of a staircase polygon scales 
with its perimeter as n. Also, it can be shown that the mean column height of a staircase 
polygon scales with its perimeter (width) as n 1 / 2 . 

Higher moments of the area random variable X n can be computed similarly to the 
discussion above. For the fc-th factorial moment of X n , where k G N, we have 

E f \ [£ n ] ~rr~fcG(t) qj 
.,. , , m { m )kPn,m q q=l fa A\ 

£ mP „,„ mui • (3 ' 4) 

where (a)* = a(a — 1) • . . . • (a — k + 1). Derivatives of G(t, q) can be obtained recursively 
from the functional equation (|3.1|) by repeated differentiation w.r.t. q. The result (see 
Lemma \A. II in the next section) is given a follows. The derivative in (J3.4J1 is algebraic, with 
leading singular behaviour 



1 d k 

k^ GM 



lk ' 1 • Oil/ r /)' -)). (3.5) 



where 7^ = 3fc/2 — 1/2, and /& (strictly) positive. A recursion for the coefficient is given 
by 



16 

0<Kfe 



where fceN and / = —1. Inserting ()3.5|) into f)3.4|) and extracting coefficients, using ()3.3|) 
and standard transfer theorems ^U], gives 

E[(jf n ) fc ] = H - ^ 7 °| / fc n 3fc / 2 (1 + O^ 1 / 2 )) . 

/oC 7 ° r (Tfc) 
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The (factorial) moments of X n diverge for large perimeter. In order to obtain a limit 
distribution, we introduce a normalised random variable X n of area by setting X n = 
X n /n 3 / 2 . This normalisation is natural in view of the scaling of the mean area with the 
perimeter, as described above. Alternatively, the area is the sum over diagonal lengths 
(column heights) and thus scales like n ■ n 1//2 , as argued above. The moments of X n are 
then asymptotically given by 

E[X n fc ] = E[(X n ) k ] (1 + 0(n-V>)) = J-^f^h (1 + 0(n-^)) , 

where we expressed the k-th moments in terms of the r-th factorial moments, r = 0, 1, . . . , k. 
The recursion ()3.6|) can be used to show that the values rrik = tim^-^ E[X*] satisfy the 
Carleman condition ^fcC 771 *) -1 ^ = +°°> implying the existence and uniqueness of a limit 
distribution with moments rrik This distribution is known as the Airy distribution, see 
[TTj for a review. For models of simply generated trees, corresponding convergence results 
appear in [TJ [T3] . 

We finally describe a convenient technique to compute the numbers fk, which will be 
central in the following sections. This approach, called the method of dominant balance 
|34j . was initially followed when analysing the scaling function of staircase polygons |57| . 
Whereas it has been applied under a certain analyticity assumption on the corresponding 
model, it can be given a rigorous meaning at a formal level, see the remark after Lemma 
14. II in the following section. It has been proved [331 Thm. 5.3] that there exists a uniform 
asymptotic expansion of the perimeter and area generating function G(t, q) about t = t c , 
which is, to leading order, given by 

G(t, q) ~ G(t c ) + (t c - t)^F ( (te 1 _~^ 3 / 2 ) (* "> *« ), (3-7) 

uniformly in q < 1 near unity. The function F (e) is called the scaling function and admits 
an asymptotic expansion about e = 0, the exponents 1/2 and 3/2 are called the critical 
exponents of the model, see also [HU UH] ■ The scaling function is explicitly given by 

K ' Ai(2 8 / 3 e- 2 / 3 ) 64 8192 262144 v ; ' 

where Ai(z) = - J °° cos(t 3 /3 + tz)dt is the Airy function. The coefficients fk of ()3.5|) 
appear in the scaling function via F (e) = J2k>o(~^) k f kek - ^ n ^ s ^ s seen ^ asymptotically 
evaluating the lhs of ()3.5j) . using the scaling form 1)3.7)1 . and by comparing to the rhs of 
(13.5)1 . The scaling function can be extracted from the functional equation (j3.1j) as follows. 
Insert the scaling form (J3.7)) for the generating function G(t, q) in the functional equation 
(13.1)1 and introduce new variables s and e by setting x = x c — s 2 and q = 1 — es 3 . Then, the 
functional equation yields a differential equation for F (e), when expanded up to order s 2 , 
see also [SH] for a detailed account. The differential equation for F (e) is of Riccati type 
and given by 
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This immediately translates into the above recursion ()3.6|) for the coefficients fk of the 
scaling function. 



4 Limit distributions and scaling functions 

We now consider the case of arbitrary M. Whereas this leads to cumbersome expressions, 
we are following the same ideas as in the previous section. Recall that the set of staircase 
polygons of half-perimeter no is denoted by V no . We consider a uniform model where 
each polygon p G V no occurs with equal probability. Then, to each counting parameter 
a corresponding random variable is attached. This leads us to consider discrete random 
variables Xk, no for the diagonal length moments rik(p), for k — 1, . . . , M. The generating 
function G{u) is a formal power series in u, i.e., G(u) G C[[tt]], where . . . , Xjv]] 

denotes the ring of formal i?-power series in the variables X\, . . . ,x^. In fact, G{u) G 
Z[ui, . . . , um] [[ m o]]? where R[xi, . . . , x^} denotes the ring of i?-polynomials in the variables 
Xi, . . . , xn- The function G(u ), where u = (u , 1, . . . , 1), as considered in the previous 
section, is the (half-) perimeter generating function of the model. Write the generating 
function G{u) as 

= Pn ,n 1 ,...,n M Uo°U 1 1 • . . . • U A f . 

n ,ni,...,n M 

The coefficient p no , ni ,...,n M is the number of staircase polygons of half-perimeter n where, 
for k = 1, . . . , M, the k-ih diagonal length moment has the value n^. We have < 
E ni ,.., nM P%,ni,.,nu = [ u o°]G{u ) = C„ _i < oo, where C n is the n-th Catalan number. 
In a uniform model, where each polygon p G V no has the same probability l/C no _i, the 
discrete random variables Xk, no corresponding to the diagonal moments nk{p) have the 
probability distribution 



P(^i,n — n-i, . . . , X Mt n — n M ) 



Pno,ni,...,n M 



Sni,...,n M Pn ,ni r 



We will analyse the joint distribution of the random variables Xk >no , where k — 1, . . . , M, 
in the limit of large perimeter. This will be achieved by studying asymptotic properties of 
mixed moments of the random variables Xk, no i given by 



V n kl ■ r) ku r) 

k X yk M i _ l^n u ...,n M U l ••• M ^0,"i "AJ 

,n ' • • • ' ^-M,n Q \ 

Z^m,...,n M Pn ,ni 



HKno ■■■■■ X 



n M 



We employ a generating function method. Note that we have 

fcl / „ \ fcjw 



^L A l,no ' • • • ' A M,n J 



K1 {u x £-y .....(u M ^) M G(u) 



U=UQ 



K°]G(u ) 



We can thus study asymptotic properties of the moments by studying asymptotic prop- 
erties of the coefficients of the generating functions appearing on the rhs of the above 
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equation. The latter problem can be treated by singularity analysis of generating func- 
tions, as indicated in the previous section. It will prove convenient to consider factorial 
moments, given by 



E[{X 1>no ) kl ■ . . . ■ {X 



M,n )k M \ 



U=UQ 



K°]G(u ) 



Moments can be expressed linearly in factorial moments. It will turn out below that 
they are, for ki, ...,kM fixed, asymptotically equal. Let us introduce factorial moment 
generating functions gk{uo), defined by 



\ Qkx Qk M 



9h(uo) ■= ^^■■■— j ^G{u) 



U=UQ 

where k = (ki, . . . , ku) £ Nq 1 . We use the multi-index notation. In particular, k\ = 
k\\ ■ . . . ■ ku 1 --, \k\ = k\ + . . . + kj\<i, and k < I if ki < li for i = 1, . . . , M. We will also use unit 
vectors e^, k = 1, . . . , M, with coordinates (ek)i = b~i,k for i = 1, . . . , M. We clearly have 
gk(u ) e C[[uq\] for all k. The function g (u ) = G(u ) = G(u c ) + Y,Zo foA u c ~ «o) (m)/2 
has been studied in the previous section (|3.2|) . where we found u c = 1/4, G{u c ) = 1/4, 
/o,o — — 1? fo,i — 1, an d /o,z = for / > 1. The functional equation (|2.2|) can be used to 
show that the properties of go(u ) carry over to those of gk(uo) f° r k ^ 0. We have the 
following lemma. 

Lemma 4.1. For k ^ 0, all factorial moment generating functions gk{uo) are algebraic. 
They are analytic for \uq\ < u c = 1/4, except at uq = u c , with Puiseux expansions of the 
form 

oo 

^K) = ^/ fc ,/K-n )'/ 2 -^, (4.1) 
z=o 

where 7& = — 1/2 + + i/2)ki. The leading coefficients fk,o = fk are given by 

fu = Ckfa fel • • • • • fell j where the numbers Ck are, for k 0, determined by the 
recursion 

M-i ^ 

C k = -27fc_eiCfc- ei + ^ {h + l)c fc _ ei+1+et - - QC fc _f, (4.2) 

i=l l^O,l^k 

0<l<k 

with boundary conditions Co = 1 and Ck = if kj < for some 1 < j ' < M. The coefficients 
fk are strictly positive for k ^ 0. Explicitly, we have f = — 1 and f £k = k\2~ 3( - k+1 ^ for 
k = l,...,M. 

Proof. The above lemma is a special case of [HJ Prop. 5.3, 5.6]. We outline the idea of 
proof. 
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Successively differentiating the functional equation ()2.2j) yields gk(v>o) as rational func- 
tion of go(uo) and its derivatives. Due to the closure properties of algebraic functions, 
these functions are algebraic again. The functions gk{uo) then inherit analytic properties 
of go(uo). In particular, gk(uo) is analytic for \uq\ < u c , except at u = u c , with a Puiseux 
expansion of the form (j4.1j) for some exponent jk- The value of the leading exponent 7*. and 
the recursion for the coefficients fk (hence for the coefficients c&) is proved by induction, 
via asymptotically analysing the fe-th derivative of the functional equation. □ 

Remarks. 

i) For the model (|2.4|) of staircase polygons, counted by width and (vertical) column height 
moments, Lemma 14. II applies as well, however with different numbers u c = (1 — y^) 2 for 
< y < 1, fo — — y -1 / 4 and f ek = k\2~ k y~ k l A for k = 1,...,M. It has been shown 
that Lemma 14. II applies in a much more general context |4"T| Prop. 5.6]. In particular, the 
amplitude ratios /fe/o fc ' 1 fe 1 kl • • • • • fe!l M = c k are model independent in that situation, i.e., 
independent of u c , f and f ek . 

ii) The recursion ()4.2|) for the coefficients fk can be obtained mechanically, if a tight bound 
on the exponents 7& is known. (Typically, analysis of the first few functions gk{uo) suggests 
the form of 7*.. Then, one can prove that 7*, is an exponent bound by induction, using 
the functional equation. Such a proof is generally easier than a full asymptotic analysis 
of the functional equation, see [JTJ Prop. 5.5].) If we replace the functions gk(u ) in the 
expansion of G(u) about u\ — . . . — Um — 1 by their Puiseux expansions ()4.1|) . we then 
have 



G(u) = ^(-l)' fc W«o)(l - «i) fcl • • • • • (1 - u M ) 



k (4.3) 
= G(u e ) + V^ToF ( l _ ~^ l+2)/2 , . . • , l _ -J^ +2)/2 , Vu-^To^ 

where, for e = (ei, . . . , ejv/), the function F(e, s) G C[[e, s]] is a formal power series, and 
F(e, s) = Fi(e)s l , where Fi(e) = J2k(~^y k ^fk,i ek - 111 statistical mechanics, the function 
F Q (e) is called the scaling function, and the functions Fi(e) are called the correction-to- 
scaling functions, if the equation ()4.3j) is valid as an asymptotic expansion in u c — Uq, 
uniformly in 1 — q\, . . . , 1 — qM, see For a probabilistic interpretation, see the remark 
after Theorem 14.21 

The functional equation (j2.2j) induces partial differential equations for the functions 
Fi(e). They are obtained by substituting ()4.3|) in ()2.2|) . and by introducing variables 
s,€i, . . . ,€m via uq = u c — s 2 and = 1 — €kS k+2 for k = 1, . . . , M. Then, a partial 
differential equation for Fi (e) appears in the expansion of the functional equation in s at 
order s l+2 . Computing these equations is called the method of dominant balance [HI]; see 
for the case M = 1, and jH] for the case of general M. For example, we get for F (e) 



the equation 
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This yields, for k ^ 0, a recursion for the coefficients fk of Fo(e) of the form 

1 M_1 j _L 1 

Y^7fc-e 1 /fc-ei + -^—{h + l)/ fe _ ei+1+ei + ^ flfk-l = 0, 

1=1 0<Z<fc 

which translates into the above recursion ()4.2|) for the coefficients Cfc. 

We are now prepared for a discussion of moment asymptotics. Recall that the random 
variables X k ^ no are related to the diagonal height moments n k (p), which are defined in 
(|2.1|) . Since the number of diagonals of a staircase polygon equals its half-perimeter minus 

1/2 

one, and the average diagonal length of a staircase polygon scales like n ' , we expect that 
the mean fc-th diagonal length moment scales like n]^ k ^ 2 . We thus introduce normalised 
random variables X k>no by 

(Y Y Y \ — I ^ l ' n ° ^ 2 '"Q X M ,n \ (A p-n 

^l,no)^-2,n ) • • • i^M,n ) I 3 / 2 > 4/2''''' (M+2)/2 I' V ' 

Then, the mixed moments of the random variables X k ^ no converge to a finite, nonzero limit 
as n -> oo. We have the following result. 

Theorem 4.2. The random variable (X l no , X 2 . no , . . . , X^no) of eqn. ()4.5|) converges in 
distribution, 

(Xi tno , X2 )no , . . . , Xm,u ) ~^ (otiXi, . . . , OLmXm ), 

where X k = J e k (t)dt denotes the integral of the k-th power of the standard Brownian 
excursion of duration 1, and the numbers a k are given by a k = — We also have 

moment convergence with 

eia- ■ . . . ■ xfcj = (i + <W: 



where the numbers u c , 7a, and fk are those of Lemma \4-1\ The numbers a k are explicitly 
given by a k = 2 _9fc//2 . 

Proof. Theorem 14.21 is obtained by specialising Thm. 8.5] and [JTJ Cor. 9.1] to the 
functional equation ()2.2|) . We sketch the proof. 

By standard transfer theorems jTHj for algebraic functions, the results of Lemma 14.11 
translate to an asymptotic expression for the factorial moments. It is 

E[(X hno ) kl ■■■■ (X M , no ) kM ] = ^.^l°\ f k n^' + 0(n^ 1/2 ), 

JqUc 1 ilk) 

The explicit error term implies that factorial moments are asymptotically equal to the 
(ordinary) moments E[X^ • . . . ■ Xj^ 1 ]. In particular, the moments EfX^ • . . . ■ Xj^ 1 ] 
converge to a finite, nonzero limit as no —>■ oo. We have 

¥lxkl yk M i _ E [ X l*n • • • • • X M,n \ _ k\ r(7p) / n(n~ 1/2 \\ 

nXl,n ■■■■■ X M,n \ ~ " /oM 7 fe -7o r(7fc) /fc [} + <?( W )) ' 
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Set rrik = hin^^oo E[X^ o • ... • X^ 1 ]. The recursion (|4.2|) or the partial differential 
equation (|4.4jl induces an upper bound on the growth of the coefficients This can be 
used to show that for t £ R 

lim — — — = U. 

|fe|-voo fc! 

This guarantees the existence and uniqueness of a limit distribution with moments rrik 
via Levy's continuity theorem j^j. The connection to Brownian excursions follows by 
comparing the moments, see |HJ Thm. 2.4]. □ 

Remarks. 

i) For a discussion of the case M = 1, see the previous section. The case M = 2 has been 
analysed in [301 EI]. The above approach yields the existence and uniqueness of a (joint) 
limit distribution for M arbitrary, and provides recurrences for their mixed moments. It 
also allows to analyse corrections to the limiting behaviour. It applies in a more general 
setup of models whose generating function satisfy a g-functional equation This in- 
cludes simply generated trees as a special case Pfl 113] . For M = 1 and M — 2, moment 
recurrences for Dyck paths have also been obtained [21] via Louchard's formula [27] from 
stochastics. The function Fo(e) in ()4.3)1 is a certain Laplace transform of the moment 
generating function for the limiting probability distribution, see [SHI Thm. 7]. 

ii) For the model ()2.4j) of staircase polygons, counted by width and vertical column height 
moments, Theorem 14.21 applies as well, however with different numbers u c , fo and f ek , see 
the remark after Lemma 14.11 The numbers are, for k = 1, . . . , M, explicitly given by 
«fc = j / (*- 1 )/ 4 2 3 - 5 */ 2 . The same remark also implies that the ratios 



HKn • • • • • ggj = fc!c fc r( 7o ) 1 -| fc ir (7ei )^.....r( 7eM )^ , , 1/2 



E[X hno f ■ ...■E[X M ,n } k ^ r( 7fc 

(4.6) 

are asymptotically model independent, i.e., asymptotically independent of u c , fo, and 
fe k , where k = 1, . . . , M. In particular, the corresponding numbers coincide for staircase 
polygons, counted by perimeter and diagonal length moments, and for staircase polygons 
counted by width and column height moments. For this reason, we will use the ratios ()4.6)1 
to check for a particular type of limit distribution in the following section. 



5 Universality and self-avoiding polygons 

For self-avoiding polygons, the notion of diagonal length moments as studied in the previous 
sections is not well-defined, since the intersection of a diagonal line with the (area of the) 
polygon need not be connected. However, we can consider a self-avoiding polygon p of finite 
perimeter as being built from a set of diagonal layers DL(p), where each layer consists of 
a finite number of connected segments d, see Figure If we restrict to staircase polygons, 
each layer consists of a single segment, see Figure |U We define parameters n£ (p) and 
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4 (p) by 

4 a) (p)= E few) > 4 6) (p)- E few*), 

L<=DL(p) \d£L / LeDL(p) \d£L / 

where is the length of a segment d of a diagonal layer L, see also Figure El Let us 
call these parameters k-th diagonal layer moments. For staircase polygons p, we have 

6 7 
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Figure 5: A self-avoiding polygon with eight diagonal layers, which each consist of one or several 
connected segments. The second and sixth diagonal layer from below each consist of two diagonal 
segments. The other diagonal layers each consist of a single segment. 

n k(p) = n i a \p) = n t (p)- For self-avoiding polygons p, we have nf\p) = nf\p), and 
nS a '{p) is the area of the polygon. The case k — 1 has been studied in a number of 
different investigations |37||3H], supporting the assumption that the limit distribution of 
staircase polygon area and self-avoiding polygon area are of the same type. Apparently, the 
above parameters have not been studied previously for k > 1. We estimated moments of 
the parameters n^\p) and nf\p) by a Monte-Carlo simulation of self-avoiding polygons, 
within a uniform model where each polygon of fixed perimeter occurs with equal probability. 
The algorithm is described in [28J. Each Monte-Carlo step consists of an inversion or 
a certain reflection of a randomly chosen part of the self-avoiding polygon. Polygons 
which are no longer self-avoiding are rejected. Checking for self-avoidance is done in time 
proportional to the polygon length. Determination of the segment structure of diagonal 
layers is done "on the fly" in time proportional to the polygon length. For a given half- 
perimeter n , we took a sample of 10 6 polygons, with at least 10 x n Monte-Carlo update 
moves between consecutive measurements. We used the random number generator ran2 
described in |2> • Estimates for the r-th moment of the random variable corresponding to 
the parameter njj. are, for fixed half-perimeter no, given by 

m k,n V r ) ■- V ~ 

where p no , n denotes the number of sampled self-avoiding polygons of half-perimeter n and 
value n of the parameter nj^ . An analogous expression is used for the r-th moment m£ „ (r) 
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of the random variable corresponding to the parameter n k . If the parameters njf , n h , and 
rtk have limit distributions of the same type, then the moment ratios TUk n {r)/m^ nQ {l) r ^ 

m k^n ( r ) / m k^n (^Y an d ^[^fc,n ]/^[^,«J r are asymptotically equal, see the remark after 
Theorem 14.21 

For comparison, we present data of the case k = 1 for the second area moment. Its 
asymptotic moment ratio is, for staircase polygons, given by 




0.0025 0.005 0.0075 0.01 0.0125 0.015 



Figure 6: Plot of the moment ratio m^ Q (2)/m^ o (l) 2 against l/(2no). The line is a least 
square fit through the data points. 

and extrapolated the asymptotic moment ratio by a least square fit, obtaining the value 
1.06084, see Figure El Both the value of the least square fit and the small data spread in 
the figure indicate that this value is consistent with the corresponding number for staircase 
polygons. 

For k = 2, we present data for the second diagonal layer moment in Figure Q using 




0.0025 0.005 0.0075 0.01 0.0125 0.015 Q 0.0025 0.005 0.0075 0.01 0.0125 0.015 



Figure 7: Plots of the moment ratios (2) /m^^il) 2 (left) and rn2.n (^)/ m 2,n (^) 2 (right) 
against l/(2no). The lines are least square fits through the data points. 
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the two different definitions n 2 and n 2 . We sampled self-avoiding polygons for perimeter 
values 2n G {64,128,256,512,1024}. Extrapolating the asymptotic moment ratios by a 
least square fit yields the values (2) / V?"4n (l) 2 — > 1.2162 and rnfl lo (2)/mfl lo (l) 2 — ► 
1.3088. For staircase polygons, the corresponding ratio is 



This indicates different limit distributions for both models, in contrast to the case k — 1. 

We also considered a corresponding generalisation of column height moments, which 
we call vertical layer moments. In section |21 we analysed the model of staircase polygons, 
counted by width and column height moments. We found that this model has a limit 
distribution of the same type than that of the model with diagonal moments. We first 
analysed the limit distribution of the column height moments for the model of staircase 
polygons counted by total perimeter. An exact analysis of the first few moment generating 
functions for M = 2 and M = 3, using the functional equation ()2.4j) . is consistent with 
the assumption that the limit distributions of both models are of the same type. We 
also performed a Monte-Carlo analysis on self-avoiding polygons for the corresponding two 
types of vertical layer moments, where we sampled polygons w.r.t. half-perimeter n Q . We 
found that the moment ratios (2) / 'm^Jl) 2 and n (2) / m 2n (^) 2 f° r vertical layer 
moments yield (within numerical accuracy) the same values than the corresponding values 
for diagonal layer moments. This suggests that the limit distributions are of the same type, 
being however different from those of staircase polygons. 



We consider some classes of directed square lattice random walks. These have been analysed 
in j^Oj by a generating function approach, mainly according to their area laws. We extend 
this analysis by providing moment recurrences for the laws of counting parameters of 
generalised area. A main ingredient in our approach is the method of dominant balance. 

All walks start at the origin. The only allowed steps are forward unit steps along the 
positive and negative diagonals. Such walks are called Bernoulli random walks. If the 
walk does not step below the horizontal axis, it is called a meander. A Dyck path is a 
meander terminating in the horizontal axis. A bilateral Dyck path is a Bernoulli random 
walk terminating in the horizontal axis. These walk models are discrete counterparts of 
Brownian motion, meanders, excursion and bridges. Corresponding convergence results 
appear e.g. in [221 El EI] ■ 

For a Bernoulli random walk b, let its length be the number of its steps n(b). We 
are interested in the k-th moments of (the absolute value of) height, defined by nk(b) = 
J2 S £w IM s )| fe > where h(s) is the height of the walk at position s, with s = 0, 1, . . . ,n(b). 
Define the weight W{,(it) of a Bernoulli random walk b by 



lim 



E[(X 2 , no ) 2 ] = 19 
E[X 2 , no P 15 



1.26667 



6 Related models 



w b (u) = u, 



n{b) 



ni(6) 



" • • • ' U M 



n M {b) 
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and let G"(u) = J2beB Wb ( u ) denote the generating function of the class B of Bernoulli 
random walks. Likewise, define generating functions for the other classes of random walks 
by restricting the summation to the corresponding subclasses of Bernoulli random walks. 

Theorem 6.1 ([30J). Let G^(u), G®(u), G {m \u) and G^(u) denote the generating 
functions of Dyck paths, bilateral Dyck paths, meanders, and Bernoulli random walks. The 
following functional equations are satisfied. 

1 

l-ul Ul -...-u M GW(v(u)y 
1 

l-2ul Ul - ...-umG^^u))' 
GW(u)(1 + m ....-w m G^(u(u)), 
G^(u)(1 + 2u - ...-umG^^u)), 

where the functions ffc(ii) are given by 

M n 

v k {u) = Y[u} k) (k = 0,l,...,M). 

l=k 

□ 

A proof of these formulae has been indicated in [SO]- The underlying combinatorial 
constructions are as follows. Dyck paths are ordered sequences of arches, where an arch 
is a Dyck path, which does not touch the horizontal line, except for its start point and 
its end point. Bilateral Dyck paths are ordered sequences of positive or negative arches. 
Meanders are either Dyck paths or Dyck paths, followed by a meander with an additional 
base layer attached. Bernoulli random walks are either bilateral random walks, or bilateral 
random walks followed by a positive or negative meander, with an additionally attached 
base layer. These constructions translate immediately into the above functional equations, 
where the same techniques as in the proof of Theorem 12.11 are used. See also jH] for Dyck 
paths. 

Singularity analysis of the moment generating functions can be done in analogy to stair- 
case polygons. Bounds on exponents can be proved by induction, using the corresponding 
functional equations. Recursions for coefficients can then be obtained by applying the 
method of dominant balance. We have the following lemma. 

Lemma 6.2. All generating functions g\)(uo) are algebraic, where (■) G {(d), (b), (m), (r)}. 
They are analytic for \uq\ < u c = 1/2, except at uq = ±u c , with Puiseux expansions about 
u = u c of the form 

oo 

9i } (u )=j2& u *- u °y /2 ^- 

1=0 



G {d) (u) = 

G^(u) = 

G {m) (u) = 
G {r) (u) = 
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The exponents 7^ are given by 



(r) (d) . <3 

"7fc 



7fc = ~2 + 2^ V 2 J *' Th =Th =Th +h 7fc 

i=i ^ ' 

TTie leading coefficients fj^ Q = fj^ are, for k ^ 0, determined by the recursions 

Af-l 

7^/^ +2^ + 1)^+ l)/^ ei+1+ . + E /^/^ = 0, 

i=l 0<i<fc 
Af-l 

7^/^ + 2 + !)(** + l)/^ ei+1+ei - 8 £ = 0, 



i=l 0<i<fc 
Af-l 

(m) ,(m) . r> \~~V„- 1 1 i\-f( m ) 1 fHfW 



7^/^ + 2 £ (i + 1)^ + i)ftl +1+ei + E /rvri = 0, 

i=l 0<i<fc 

,(r) _ f (6) f (m) 

Jk ~ 2—t Jl J k-li 
0<l<k 

with boundary conditions /g = —4, /<j = 1/2, /q'™' 1 = 1, /g r ' ) = 1/2, and /£' = if 
kj < /or some 1 < j < M. The coefficients f^} are strictly positive for k ^ 0. □ 

Remark. The relations between the coefficients for Dyck paths, bilateral Dyck paths, 
meanders and Bernoulli random walks have discussed in [HO] in the case M — 1, and in 
[3T] partly in the case M — 2. The recurrence for Dyck paths for general M is given in 

SB- 

The above recursions can also be phrased in terms of the corresponding coefficient gen- 
erating functions Fq (e) = J2k(~ l)' fe '/fc^ efe 5 which are (candidates for) scaling functions. 
They are given by 

(, Ms . x \ Af-l 

~F W (e) + + 2 E( j + W^(«0 + = 16, 

i=l ^ ' ' / i=l 

F^(e) = -2 (Fo^Ce))" 1 , 

. (i^(e) + g (l + I) ,^-(e)) + 2 g (* + («) + 

+ F (m) (e)F W ( e )=4, 
F (r) (6)=F (6) ( e )F (m) (e). 



(6.1) 



As in section one can now proceed in defining random variables and prove, after 
appropriate normalisation, the existence and uniqueness of a limit distribution. For Dyck 
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paths and k-th height moments, the limit distributions coincide with those of k-th excur- 
sion moments, i.e., the integral J e h (t) dt over the k-th power of the standard Brownian 
excursion of duration 1. In view of the bijection described in the introduction, this result 
is not unexpected. For meanders, bilateral Dyck paths and Bernoulli random walks, we 
have convergence to k-th moments of meanders, of the absolute value of Brownian bridges 
and Brownian motion. In fact, convergence results for underlying stochastic processes have 
been obtained previously [221 El EH- Note that, apart from moment convergence, the above 
method also yields moment recurrences and corrections to the asymptotic behaviour, which 
cannot easily be obtained following the stochastic approach, see jH] for the case of Dyck 
paths. Consistency with the stochastic description has been demonstrated for k = 1, by 
deriving explicit expressions for the scaling functions 30J. Also, the second and the last 
equation of (j6.1J) have a stochastic counterpart, see [301 Thm. 2] and [301 Thm. 7]. We 
leave further details to the interested reader. 

7 Conclusion 

We analysed the model of staircase polygons, counted by perimeter and k-th diagonal 
length moments, and by width and A;-th column height moments. The model may serve 
as a toy model of vesicle collapse in k + 1 dimensions. For example if k — 2, staircase 
polygons may be converted into a three-dimensional object (vesicle) by replacing each 
column of height n by an n x n slice extending into the third dimension. The second 
column height moment is then the volume of the vesicle. 

We derived limit distributions for the diagonal length moments in the limit of large 
perimeter. These results may be viewed as a special case of limit distributions correspond- 
ing to g-functional equations with a square-root singularity as dominant singularity of the 
perimeter generating function jUj. Distributions of the same type appear for the column 
height moments in the limit of large width. They also seem to appear for column height 
moments in the limit of large perimeter, but a proof is an open problem. 

We presented corresponding results for related, but different problems of directed walks 
such as meanders and bridges. Variants of these models, together with a contact activity, 
may be useful for questions of polymer adsorption [T^]. Also, the methods used in this 
paper can be applied in the analysis of parameters on trees, such as the Wiener index [20] 
and left and right pathlengths [21], see also [H]. 

It is perhaps surprising that staircase polygons and self-avoiding polygons might obey 
the same area law, but that their higher layer moments seem to be distributed differently. 
Apparently, both models do not share the same universality class. Understanding self- 
avoiding polygons remains a challenge for the future. 
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